Nucleon to Delta Weak Excitation Amplitudes in the Non-relativistic Quark Model 



Jim Liu, Nimai C. Mukhopadhyay and Lisheng Zhang 
Physics Department, Rensselaer Polytechnic Institute, Troy, NY 12180-3590 

We investigate the nucleon to Delta(1232) vector and axial vector amplitudes in the non- 
relativistic quark model of the Isgur-Karl variety. A particular interest is to investigate the SU(6) 
symmetry breaking, due to color hyperfine interaction. We compare the theoretical estimates to 
recent experimental investigation of the Adler amplitudes by neutrino scattering. 
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I. INTRODUCTION 

Recent claim by a group of experimentalists at Los Alamos that they may be observing neutrino flavor oscillations 
is a dramatic example of high topical interest in the low and medium energy neutrino physics. Such experiments are 
now possible, for example, at Brookhaven, Fermilab and Los Alamos, just to name a few laboratories among many 
facilities around the world. As a by-product of this experimental opportunity of having excellent medium energy 
neutrino beams [E v ^1 — 3 GeV) being readily available, a rebirth of exploring with neutrinos the properties of 
hadrons, particularly nucleons, both in their ground states and the resonance region, is expected. There is a long 
history of such investigations [2-18] in the Delta(1232) resonance region since the sixties, the most recent experiment 
p8| in the Delta region have been done at Brookhaven. Given the topical interest of the structure of hadrons from 
£vq the QCD point of view, the exploration of nucleons and their excited states, by both electromagnetic and weak probes, 
merit special attention. In this exploration, the weak structure functions are difficult to determine, but they provide 
t-H valuable information, often complementary to that obtained by the electromagnetic interaction. 

Among all the excited states of the nucleon, the Delta(1232) is perhaps the best studied one ]l9| , p0| |, by strong, 
electromagnetic and even weak interactions, the last one being the focus of this paper. Along with the nucleon, the 
Delta is of fundamental interest to quantum chromodynamics (QCD) and its application to the problem of hadron 
structure. Until we learn how to use QCD rigorously to compute hadron properties, particularly in its low-energy 
(hence non-perturbative) domain, we must cope with models that are "QCD-inspired" . One of the most successful 
I/-) , examples of these imperfect constructs is the quark shell model (QM), in which the gluonic degrees of freedoms are 
replaced by effective potentials. The origin of this can be traced back to the pre-QCD sixties, when Giirsey, Radicati 
p0[ and Sakita |2l[] described the nucleon and the Delta in the fundamental 56-dimensional representation of the 
spin-flavor SU(6) symmetry group. In this limit, the nucleon and the Delta wave functions are \N >— \{56} 2 S S > 
and |A >= {56} 4 S S > . They are degenerate in the symmetry limit. The degeneracy is lifted by the color hyperfine 
^H. interaction [22|. Thus, in the Isgur-Karl QM, the wave functions of the nucleon and the Delta are [^3|,|4) : 

\N >= a s \ 2 S s > +a s ,\ 2 S s , > +a M \ 2 S M > +a D \ i D M > +a p \ 2 P A >, (1) 
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|A >= 6 S | 4 5 S > +b s ,\*S s < > +b D \ l D s > +b D ,\ z D M >, (2) 

where the a's and b's are determined by diagonalizing the QM Hamiltonian in the N=2 harmonic oscillator basis. 
The purpose of this work is to investigate the process |2) 

p+v^ A++ + /X", (3) 

in the framework of the non-relativistic QM, using wave functions (2). Despite a long history of theoretical (Refs.[2- 
14]) and experimental [15-18] investigations of the process, no calculations of the amplitudes for the weak transitions 
are available in the literature in the context of these general wave-functions. Our main objective in this paper is to 
remedy this. We shall focus in this paper on the all relevant helicity amplitudes. There are four transverse ones: 
A\, 2 , Ay 2 and Ay 2 , A^ 2 , where V stands for vector and A represents axial vector. We shall also discuss the relevant 
longitudinal (scalar) weak amplitudes. From the electromagnetic processes p5[ 

N + 7 t± A, (4) 

we know the vector helicity amplitudes , both extracted from the experiments and as computed in various theoreti- 



cal approaches, such as the quark model (QM), both non-relativistic [|23 24 1 and relativized [|26|, bags |27|], topological 
p3| and non-topological |2l| solitons, and by the lattice gauge theoretic method p9fl . Similarly thorough theoretical 
investigations on the axial vector amplitudes have been suggested in a recent investigation p4| in the framework of 
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QM. This is our objective here: to compute the vector and axial vector amplitudes in the framework o f QM . We shall 
relate them to the recently extracted Adler 0] amplitudes from the Brookhaven neutrino experiment Jl§| | , which has 
studied the reaction 

v^ + d^fx- + A++(1232) + n s , (5) 

n s being spectator neutron. When necessary, we shall go back to the older experiments as well. 

Our main goal is to investigate the success of the non-relativistic QM in reproducing the measured transition form 
factors for the process (3). In the case of (4), the phenomenologically extracted magnetic dipole (Ml) and electric 



quadrupole (E2) amplitudes |25| at the real photon point are considerably larger than the values obtained in the 
QM. Their values away from the real photon point are not very well-known experimentally as yet. This is going to 
change with the advent of the electron facility called the CEBAF. The theoretical deficit of the transition magnetic 
amplitude, compared to the observed one, seems to be confirmed by the low-energy Compton scattering, where the 
resonant magnetic polarizability appears to be completely dominated p0| by the Delta contribution. One of our 
goals here is to see if the axial-vector analogue of the Ml amplitude in the process (3) can be reproduced in the 
non-relativistic QM. We shall also make an estimate of the axial-vector analogue of the E2 amplitudes, which would 
be zero in the SU(6) symmetry limit. Thus, its non-zero value would be a direct manifestation of the effect of the 
color magnetism, just as a non-zero value of the E2 amplitude is in the vector sector. 

Remainder of this paper is organized as follows: We derive an effective Hamiltonian for the QM with vector and axial 
vector interactions in section II, wherein we also find a relation between transverse vector and axial vector amplitudes. 
The calculation of helicity amplitudes is done in section III. We find relations between the helicity amplitudes and 
Adler's form-factors in section IV, for the comparison of our model estimates with experiments. We collect some 
exact SU(6) relations in section V. Section VI contains the main results of our QM calculation and its comparison 
with experiments. In this section, we introduce an axial vector analogue of the E2/M1 ratio. The SU(6) breaking 
results in this ratio becoming non-zero. Neutrino experiments can give us an estimate for this ratio. Our conclusions 
are summarized in section VII. 

II. THE QUARK TRANSITION HAMILTONIAN 

The electroweak interaction at the quark level is incorporated in the Dirac equation in the usual fashion. We 
introduce the leptonic current for the weak process, which is the analogue of the electromagnetic vector potential 
A,: 

3 

fewV ~ m i ~ ^ (i " 75) W) « = o, (6) 

1=1 

so that the Hamiltonian is 

3 

H = (« ' (Pi - (1 - 75) 1.) + /3m, + (1 - 75) l ,i) . (7) 

i=l 

By doing a free Foldy-Wouthuysen reduction [plf, the quark Hamiltonian can be truncated to 
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1=1 



with 

Hit = « y E (*<k " T- fa + P*) - T- a ■ (V X I,) + O (m- 2 )) , (9) 



i=l 



Hit = a A J2 (' a ■ li + ^ a ■ fa fa + faP*) + K~ 2 )) 

2—1 ^ 1 ' 



(10) 



Here i is the quark index, a v and a A are factors, which, in general, can be different. Unless otherwise stated, we shall 
take these factors to be unity. We take quark masses mi to be the same, m, for up and down quarks. We can see that 
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the HY nt is the same as interaction Hamiltonian of the electromagnetic interaction except for the obvious difference 
in coupling constants. The lepton current can be defined, in analogy to the electromagnetic vector potential A^, 
as follows: 

1 = \/^E ( e ^ ik r + , (11) 



fe, A 



where we introduce a "weak" fine structure constant, aw, for emitting or absorbing a weak quantum, W, with 
four-momentum (kg; 0, 0, k). Here ito is the energy transfer for q^q^ — 0, given by 

M 2 - M 2 

K ° = ~ 255.8MeK (13) 

This is the analogue of the real photon point. We use Kq in the normalization factor under the radical sign in 



Eq.(ll,12), a practice adopted by many [£6 32|. This factor cancels out in the form-factor expression, having no 
influence on the form-factors we shall calculate. For the transverse amplitude, we use the boson polarization vector 
to be e = — (l,i,0). The vector transverse Hamiltonian can be separated into two pieces: 



H Ynt. trans. = H Y + H II> ( 14 ) 



with 




taking § = s. 

The correspondence between the usual photon case and the virtual W-boson exchange discussed above, needed for 
the process (3), can be seen by the appropriate identification of the photon variables to the W-boson variables. For 
transverse axial transition, Iq does not contribute, and the other term is 




Thus we have an important relation 

n trans. ~ ~gV 1 ' ^ ' 

where e A /e v , denote the possible coupling difference between vector and axial vector parts. If we compare directly 
with photon transitions, e v is simply e. 



III. CALCULATION OF THE HELICITY AMPLITUDES 



The calculation of the longitudinal and scalar amplitudes requires a discussion of the possible lack of current 
conservation for chiral vector and axial vector currents in the QM space chosen.. Though discussions on this issue 
have been made from time to time in the literature p,p2[, this issue is not settled in our opinion. In this paper, we 
simply avoid tackling this issue. We begin here with a discussion of transverse helicity amplitudes. The leading orders 
of these amplitudes do not suffer from the uncertainties of the lack of current conservation. 
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The calculation of the transverse helicity amplitudes for the vector current is standard 24 . 3^] . The wave functions 
of A and N, taking the effects of the color-hyperfine interaction, can be expressed in terms of SU(6) basis functions 



|A >= 0.97| 4 S. > +0.20\ 4 S' > -0.097| 4 L> S > +0.065| 2 £> M > 



(19) 



\N >= 0.95\ 2 S S > -0.24| 2 ^ > -0.201^ > -0.042| 4 L> M >, 



(20) 



ignoring the tiny contribution from \ 2 Pa > in (1). These two wave functions correspond p4| to energy levels — 
l230MeV, and En — 940MeV. The transverse helicity amplitudes for the electromagnetic current are defined in the 
photon-nucleon CM frame |24[]: 



AX 



--< A, Mj = \\HY t \"/N, A = A 7 - X N , k >, 



(21) 



where the helicity A is 3/2 or 1/2. These equations are straightforwardly generalizable to the case of our interest, 
introducing virtual W boson mediating the charged weak current that produces the reaction (3). From these quantities, 
we can define the familiar 23J24J transverse electromagnetic multipoles: 



Ml = %= ( 3A V 

2V3 V 



3/2 



/2 



(22) 



/2 



(23) 



We can also define the longitudinal and scalar amplitudes L v and S v , for which the virtual boson or photon helicity 
is zero. Clearly, these two amplitudes are related |32[: 



(24) 



a relation which can |32| be violated due to truncation of the model space in the QM. The relation (24) is a consequence 
of the conserved vector current (CVC), or equivalently, the gauge invariance of the vector current. 

The calculation of the transverse, longitudinal (and scalar) amplitudes for the vector current have been done by 
many authors [^3 24 3^]. Thus we do not discuss them here. We shall give here a brief discussion for the axial vector 
amplitudes. The longitudinal axial transition operator is given by 



! = 1 



where we are considering the axial vector term only, Eq.(10). Following (15), we have 

f 2. 



H 



AL 



2Kn 



-3a 



(3) 



cxp ( -k\ z 



(25) 



(26) 



From Eq.(10), the scalar term is 



3 1 



(27) 



Following (16), we have 



(^73 {P^ a 



(3) 



ttAS 
11 int 



J 

2K„ 



(3) 



72Vi 3) ) exp {-iyfikX^ 



(28) 



where k = k/ano, Q-ho being the harmonic oscillator parameter of the IK model, which has a value p3[ of 320MeV 
from the fitting of the nucleon spectra. Here tia and tin are the principal quantum numbers of A and nucleon. The 
longitudinal and scalar amplitudes are defined respectively as 
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L A =<A-J z = -\H AL t \N-J z = ->, (29) 

S A = - < A; J z = \\H A *\N; J z = \ > . (30) 

These amplitudes are shown in Table Q. 

We shall now use the computed values of A^ 2 an d V%A^ 2 , obtained with Hj of Eq.(14), and apply the relation 
(17) to compute the A A , 2 and y/3A A , 2 . The results are shown in Table ||. In analogy to Ml and E2, discussed earlier 
(Eqs.(22,23)), we can also compute the amplitudes (M1) A and (E2) A , in the form: 

(M1) A = A(k)f A (k), (E2) A = A(k)f A (k), (31) 

Here 



and f^j(k), /jj(fc) are model-dependent, to be explicitly given later. 

QM calculations of A^ 2 , A^ 2 , L v and S v are standard in the literature. For completeness, we have collected these 



results in our notation in Table Bourdeau and Mukhopadhyay |32| have discussed the violation of (24) in the IK 
and other quark models. Numerical results of the form-factors from IK model are shown in section VI, subsection D. 
Ratios of helicity amplitudes are shown in Figs. 1 and 2. Comparison with experiments is done through the Adler 
form-factors discussed below. 



IV. THE ADLER- RARITA-SCHWINGER FORMALISM 



The standard method in the theoretical treatment of weak interaction process in (3) follows the Rarita-Schwinger 
p3"| , j34f formalism and the N — > A transition form factors are introduced following a notation due to Adler Q and 
Llewellyn Smith [Q. Thus, the invariant matrix element is 

M =< (i-A ++ \H int \vp >= GFC ^ B j a < A++\V a - A a \p > . (33) 

v2 

where Gp the Fermi constant, and 8 is the Cabibbo angle, 

r=^7 Q (l-7 5 H (34) 
is the weak lepton current. We can decompose the invariant matrix element as p|,ppC|] 



+ + %( p *)>) fA " + c ti" + §ki)V»)», (35) 

where 

F Xa = q X J a - q a j\ (36) 

Pa is the Delta four-momentum, is the Delta vector spinor and u, the proton spinor. The C^' A are the so-called 
Adler form- factors Q. These form factors can be related to the helicity amplitudes through calculation of M projected 
in different polarizations. The factor -\/3 in Eq.(35) above comes from the isospin relations, such as that between the 
weak and electromagnetic matrix elements: 

<A++|V>>= v / 3<A+|V£ M b>. (37) 
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After separating j a , we get 



-L<A ++ |-A>>=^(f 7A + ^(P A)A )A 



(38) 



Recent experiments on the process (3) at the Brookhaven National Laboratory |L8| have been analyzed in terms of 
the Adlcr form factors. Thus, we must find relations between the helicity amplitudes, discussed above, and the Adler 
amplitudes. 

In the A rest frame, which is also the WN (or jN) CM frame, 



(Pa)o = Ma = E N + k = \/M 2 + k 2 + fc , 



(39) 



Pa = and q q\ = fcg — k . We also define the Delta vector spinors as 



^3/2 = e +X4 



(40) 



roX+ . 



(41) 



Here the polarization vectors are 



(42) 



Thus, the vector spinor 



V>3/2 = e +X+ 



V2 \ 



V2 \ 



(43) 



and so on. 

The nucleon spinors are 



.N 



E N +M& 



X ± 

g ' k y± 

E N +M 



(44) 



Thus we can derive the following relations: 



J] u 3 "1 I72 - u 4 "1 a72 T <->i 



M 2 ^4 



M 2 



(45) 



/2 



47rav 



jT7 ^3 "1 I72 - °4 i 572 <-^K "T- <-•« 



(46) 



47ra 



M° 3 + ^2-° 4 



~M 2 ~ 5 + 6 



(47) 



47TQ! 



If 



M 



M 2 



-CI 



kE N n y 

M 2 ° 5 



(48) 



G 



.4 



3/2 



Ana 



2Kn 



k + 



E N +M I M 



1 



-,ci 



koMA rA 
~M 2 ~ Li 



Ct 



(49) 



VSA A /2 = 



4ira 



w 



2K 



k2 V c , 

Epf+M J M 3 



.1 , fc M A ^ A 



(50) 



_t a ij^aw_ 



M° 3 + 



k M A 
M 2 



c? + c A -^ 



V 2 

M2° 6 



(51) 



5- 



2 V 2X 

From the relations for the vector amplitudes, we get 



M 



kM A „ A k k ^ A 



M 2 



-ct 



M 2 ° 6 



(52) 



^3/2 



2K Q M 
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by far the ?nosi important vector contribution in the N — > A transition. We also get 

k 3 



1 3/2 



/2 



2K M 2 {E N + M) 



f-iV 
°5 ) 



(53) 



(54) 



(kL v - k Q S l 



4ira 



w 



2K 



(E N - M) CX 



(55) 



By CVC or gauge invariance, the left-hand side is zero identically. Thus we get a CVC relation, C§ = 0, a relation 
which will be mildly violated by the IK quark model calculation, due to its inability to make the left-hand side of (55) 
vanish identically [p2| . 

Rearranging the axial-vector transition helicity amplitudes, we have 



4ira w ^ A 



2K 



M / 



2 (E N - M) 



(V3Af /2 -A 



A 

3/2 



We can define relations analogous to (22) and (23): 



(Ml) A ^-^={3A^ /2 + V3Af 



/2 



(56) 



(57) 



We also have the relations 



(E2Y 



2V3 



A A - VSA? 



1 3/2 



/2 



(58) 



2Kn 



C A = -V3 



M 



2 (E N - M) 



[E2) A , 



(59) 



^a w f c A + M Ak c a 



2K Q 



At 



A I aA 
/2 + ^3/2 



k I A-koiw A 



MV 2K 



3 ) 



(60) 



(61) 



The amplitudes S A and L A are connected by the identity 
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kL J 



koS Jl = -\ -k 



Ana 



IT 



2Kn 



k 2 - V 2 

~M 2 ~ Lj& 



(62) 



We get the conservation of the axial current (CAC) in the chiral limit (m n — * 0) in which this identity reduces to 



kV 



k Q S J 



0. 



We must examine, if the chiral symmetry breaking and the truncation of the model space in the QM results in the 
violation of the last relation. The answer we find below is in the affirmative. 
Using (49-52) we can derive the following relations: 



Ca 



2K 



M 2 

k 2 V 4:Tra w 



1 / 



(V3A 



A 

3/2 



(63) 



( 



2K 
2 V -liraw 



L A _ to sA 



k 2 



V r A 

M 2 ° 6 ' 



(64) 



M 2 
kM A 



2Kn 



2 V inaw 



■s- 



M 



k ok nA 
M 2 ° 6 



(65) 



Thus, Eqs. (56,63-65) complete the expressions for the four Adler axial nuclear to Delta form factors, in terms of the 
calculated helicity amplitudes. 

The partial conservation of axial current (PCAC) relation can be expressed in another way, as discussed earlier by 
Schreiner and von Hippcl [Q. From the pion pole dominance of the divergence of the axial current, taken between 
the nucleon and the Delta states, we get the induced pseudoscalar term given by the pion pole, exactly parallel to the 
weak current between nucleon states, wherein the pion pole term also yields the induced pseudoscalar term p0| , p^[ . 
Thus, 



cl 

M 2 



2V3M (ml - q 2 



at the pion pole. Here <?a is the A" 1 



p7r + coupling constant, recently redetermined by Davidson et al. |2£ 
g A = 28.6 ±0.3 . 



(66) 



(67) 



Actually, the determination of g A depends on which method we use in the analysis, thereby yielding a much bigger 
theoretical error than that is given in (67). In Eq.(66), is the pion decay constant 



f n « 0.97m T 

Taking the limit of the divergence of the axial current as m 2 
Treiman relation 



c 5 A (o) 



9aU 
2V3M' 



(68) 

and q 2 — > 0, we get the off-diagonal Goldberger- 



(69) 



This off-diagonal cousin of the well-known diagonal Goldberger- Treiman relation has not been given much attention 
in the literature, except in the recent work by Hemmert et al. Using (69), the estimate of C^(0) is 



C£(0) w 1.2. 



(70) 



The relations (62') and (69) should be simultaneously satisfied if PCAC is to be valid. 

We shall further discuss the subject of PCAC and its validity in the QM in subsection VLB. Form-factors calculated 
here are shown in Figs. 3-10. 
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V. THE EXACT SU(6) SYMMETRY RELATIONS 

Let us first start with the vector amplitudes. In the SU(6) symmetry limit [^5| , 

Ml ^ 0, E2 = 0. (71) 

The Fermi- Watson theorem implies these multipoles to be purely imaginary on top of the Delta resonance. Thus, 

Im(E2)/Im(Ml) = EMR = 0. (72) 

In the QM, all amplitudes are purely real, hence E2/M1=0 in the SU(6) limit, giving 

AV /2 (SU(6)) = V3AY /2 (SU(6)). (73) 

Using this relation in Eq.(53), we get, in the SU(6) symmetry limit, the largest vector form factor, Ccf, given by 
(Table H 



c i < st " 6 » = - Wws) AUsum - fa*"- (74) 

Also, the left-hand side of (54) vanishes in the SU(6) limit, giving 

C%(SU(6)) = 0. (75) 

Likewise, we can get 

CX(SU(6)) = -^-CY(SU(6)). (76) 

Finally, the longitudinal vector response vanishes in the SU(6) limit: 

L V (SU(6)) = S V (SU(6)) = 0. (77) 

Thus, the only non-vanishing multipole vector amplitude in the SU(6) limit is the magnetic dipole amplitude, given 
by 

2k — 

Ml(SU(6)) = -2AX /2 {SU{&)) = —A{k). (78) 
' dm 

We can similarly discuss the axial vector amplitude in the SU(6) limit. The only non-vanishing amplitude is 

(Ml) A (SU(6)) = ^A(k), (79) 



and 



(E2) A (517(6)) = 0. (80) 



The axial vector transverse helicity amplitudes are: 



±A A /2 (SU(6)) = Af /2 (SU(6)) = - 2 -A{k). 



The scalar and longitudinal axial vector amplitudes are: 



S A (SU(6)) = ^-A(k), (81) 
9 m 



L A (SU(6)) = -^A(k). (82) 



The Adler form factors become 



9 



C 3 A (Sl/(6)) = C£(SU(6)) = 0, (83) 



C£(SU(6)) =(4= + -A= — ) e" r/6 , (84) 



V3 3\/3 m 



1 />f2 _ 2 

C?(Sff(6)) = - — -— e - k l\ (85) 
3V3 Maw 



In Figs.l and 2, we plot the ratios 



r 



A b 

b _ 3/2 (86) 



/2 



as functions of Q 2 — —q 2 , b=V, A, for the IK wave functions of the nucleon and the Delta, and the wave-functions 
p2| inspired by a model of Vento et al. (VBJ) In the SU(6) limit, this ratio should be unity for both V and A 
currents. The deviation from unity is due to the SU(6) breaking interactions. In the IK case, that is from the color 
hyperfine interaction. 

The SU(6) relations described above are only approximate, since the SU(6) symmetry is broken by the color hyperfine 
interaction. Thus, the IK model wave functions and the actual experiments would violate these relations. The degree 



of these violations is an interesting question and lots of theoretical 1 25 |3G|.|37j] and experimental |38| attention are being 
given to it, since the finding of Davidson et al. |25|] that the EMR is not zero from the existing old electromagnetic 
data. We hope our study here of the axial vector amplitudes would trigger similar interest in the weak sector. 

We note that the SU(6) symmetry limit to the nucleon and the Delta wave function provides consistency with the 
requirements of CVC. Thus, the CVC requirement of current conservation is trivially satisfied. The vanishing of Cg , 
required by the CVC limit, is again trivially true. 



The CAC relation (62') is not satisfied in the SU(6) limit. However, we get good agreement with the off-diagonal 
Goldberger-Treiman relation, Eq.(69). Thus, going to the "real photon" point k = ko, we can evaluate the C^(Q 2 = 0) 
by using (84). We get 

C 5 A (Q 2 = 0)(SU(6)) = 1.17, (87) 

compared with the off-diagonal Goldberger-Treiman value of 1.2, a good agreement. We can also see that the identity 
(62) is satisfied in the SU(6) limit of QM, as it must. However, both sides are quite large at the real photon point: 



instead of the PCAC expectation of both-sides of (62) vanishing at the chiral limit, Eq.(62'). 



VI. MAIN RESULTS OF THE QUARK MODEL AND COMPARISON WITH EXPERIMENTAL 

RESULTS 

A. Comparison with other models 

An extensive review of the old QM calculations by many authors has been done by Schreiner and von Hippel [jl0| , 
who have provided a detailed comparison between theory and experiment, as available till 1973. The readers are 



referred to their paper for a discussion. In Table [V, we summarize the predictions of various form factors, vector 



and axial vector, in different models, and compare them with the new QM results reported here, as well as the recent 



estimates of Hemmert et al. 14 



B. CVC and PCAC 



We have seen earlier that the quantity (kL v — koS v ), the left-hand side of Eq. (55), vanishes in the SU(6) limit, 
thereby satisfying the CVC. Other quark model wave functions for the nucleon and the Delta, such as those of 
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the Isgur-Karl model and the VBJ model |36|], discussed below, do not satisfy this CVC constraint |3^]. The CVC 
violations in the IK and VBJ models are demonstrated in Figs. 5 and 6 where the C\ and the C$ are plotted as 
functions of -q 2 for W=1230MeV and W=1234MeV. 

As noted earlier, the PCAC pion pole estimate of C^q 2 ), Eq.(66), is not fulfilled in the QM. The IK quark 
model, for example, produces a value of C§(q 2 ) much smaller than what Eq.(66) suggests. This failure is readily 
understandable, because the inclusion of the pion-pole term in a QM that involves only quark degrees of freedom 
and no mesons, is not legitimate. However, what is surprising is that the off-diagonal Goldberger-Treiman relation 
(Eq.(69)) is well-satisfied in the SU(6) limit and by the IK quark model wave functions. Thus, in the SU(6) limit, we 
have Eq.(87); in the IK model, we get 

Ci{Q){IK) = 1.16. (89) 
This, however, occurs due to an accident. In the chiral limit 

( L A _ ko sA \ ^ Q) (90) 

V / chiral 

in Eq.(64), and C$ would pick up a PCAC contribution, not present in our QM calculations. The latter does not 
satisfy the chiral axial current conservation and does not pick up the PCAC contribution in Cg . However, these two 
theoretical inaccuracies somehow add up to the PCAC estimate in our SU(6) and IK quark model wave functions. 
Further theoretical work is needed to illuminate the nature of this happy accident. 

While the degree of violation of CVC and PCAC should be as small as possible in the quark model of excellent 
quality, this does not mean we prefer the SU(6) limit to the IK model. That is because the color hyperfine interaction 
and other dynamical considerations make the SU(6) limit inaccurate for baryon spectroscopy. We must, therefore, 
search for a model beyond that limit. Between the IK model and the VBJ models, the former is clearly superior, as 
the violation of CVC is less in the former than in the latter, and the IK model account for the baryon spectroscopy, 
while the VBJ model is not so ambitious. 



C. SU(6) breaking 



We can go back to Figs. 1 and 2 to demonstrate the SU(6) breakings in the IK J23 24 and VBJ |36| models. These 
effects are significant, though difficult to measure experimentally. In Figs. 3 and 9, the SU(6) limits coincide with 
the theoretical predictions of CVC and off-diagonal Goldberger-Treiman relation. Here, more realistic quark models 
produce mild violation of these important constraints. Such violation can arise from the model truncation effects. 
Some authors have tried to remedy these with the introduction of the form factors at the quark level [p|JTl|] . We shall 
not do that, as we believe this is not a satisfactory remedy in the spirit of QCD. 



D. Comparison with experimental results 

The experiments on neutrino scattering in the Delta resonance region are very difficult, nevertheless they have 
been done in a number of labs, ANL, CERN, Fermilab and BNL [15-18]. These experiments deal with single pion 
production in the charged current reaction in hydrogen and deuterium. In view of the extensive literature on the 
older experiments, we shall focus here on the most recent one. 

Latest experimental results |1| come from the Brookhaven National Laboratory(BNL). The analysis by Kitagaki 
et al. |L8). makes some strong assumptions, including the neglect of the background contributions and use of the 
polynomial forms for the transition form factors, prescribed by Adler [Q. The dipole form they use is 

F V (Q 2 ) = \(Q 2 )/(1 + Q 2 /M 2 ) 2 (91) 

where X(Q 2 ) is a function introduced by Olsson et al. p3|: 

\ 2 {Q 2 ) = 1 - (0.053 + 0.017Q) sin ( j^^ j ■ (92) 

Kitagaki et al. [jl8| find, from a fit to their new data, 

M v = 0.8918$ GeV - ( 93 ) 
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We demonstrate the results from Kitagaki et al. by displaying C 3 {Q 2 ), using Adler parameters M, in Fig. 3 as a 
function of Q 2 . We also plot the earlier experimental fit by Dufner and Tsai jlO), who got a good fit to the data on 
the electroproduction of pions in the Delta region, by using the form 

lOTl 2 = (2.05) 2 (l + 9(Q 2 ) 1 / 2 ) exp ^-6.3(Q 2 )^ 3 ) . (94) 

We shall also use the simpler dipole form j[o|, which fits the above data just as well: 

c ="= 2 '° 5 ( i+ ;sr- < 95) 

In the formulae above, Q 2 is in GeV 2 . We notice that the IK model does a reasonable job of describing the 

experimental data on C 3 , except at Q 2 = 0, where it falls short (Fig. 3) in comparison to experiment. This is not 
surprising, as we have encountered this deficit already in the photoproduction of pions pq| . 

We now come to the axial-vector form factors. Kitagaki et al. |18|] assume the Adler form for the dependence on 
Q 2 , and then try to fit their new data to a range of Ma- Thus, they use 



(i + q 2 /m 2 ) 2 ' 1 = 3 ' 4 ' 5 ' (96) 



where a's b's and c's are all model-dependent parameters determined in the Adler model [Q: 

c 3 (0) =a 3 = b 3 = 0, 
c 4 (0) = -0.3, 

cb(0) = 1.2, (97) 

04 = O-h = —1.21, 

b 4 = b 5 = 2. 

By fitting My and Ma simultaneously to their data, Kitagaki et al. get, with My given Eq.(93): 

M A = 0.97±g:i^ GeV. (98) 

Their results are found to be consistent with other earlier experiments and analysis. We thus compare our QM results 
with the above parameterization, where the experimental information controls the value of Ma- Our agreement with 
the empirical forms needed by the data of Kitagaki et al. is satisfactory in the IK model (Figs. 3, 9). 

Given our theoretical difficulties in computing the longitudinal structure functions in the QM, we shall compare 
our calculations of the axial- vector helicity amplitudes, A^ 2 an( i ^1/2 > with the experiments via the form factors 
C3, C4 4 , C^. In Figs. 7-9, we display this comparison. 

Some cautionary remarks are in order regarding the quality of comparison of the IK (or any other) quark model 
and the "experimental values" of the Adler form factors Q that we have displayed in Figs. 3-10. The analysis of 
the latest experimental results on neutrino scattering by Kitagaki et al. [jl8| does not separate the individual form 
factors directly by experiment. Also, it makes strong assumptions on the Q 2 dependence of the nucleon the Delta 
form factors, basically following the polynomial forms prescribed by Adler |4|. Finally, it neglects the non-resonant 
Born contributions p5| ], which may be small, but uninvestigated at this time. 

Below we give numerical expressions of the helicity amplitudes and the Adler amplitudes in the IK model. Since 
the differences between AY',^ and \f3AYjt are small, we give Ml and (M1) A instead: 



x 3/2 v °^3/2 

Ml 



71 2m" 4 ^' ) ( 2 ' 0237 + °- 0623fc2 - 0.00133I 4 ) , (99) 

E2 = — \= — A(k) ( 6.26 + 0.627A? 2 + 0.252^ x 10~ 3 , (100) 
V3 2m V ) 

,y 1 -2 / „ . „„t-2 



S v = -— =A(k)k ( 0.0996 + 0.00133fc ) , (101) 
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L v = V3A(k)— (-0.00104 + 0.000411fc 2N ) , (102) 

771 V / 



\A 1 / . „ „„„_-2 „-r4 



(Mir = ~^=M k ) (2.0204 + 0.0625/c - 0.00133fc ) , (103) 
V3 



(^2) A = ^ A (k) (0.325I 2 - 0.251I 4 ) x 10~ 3 , (104) 

L A = -A(fc) (o.8235 + 0.0225fc 2 - 0.0009351*) , (105) 

S A = Jlh-AQi) (o.338 + 0.00546I 2 - 0.000205I 4 ) . (106) 

The IK model parameters have some uncertainties due to variable estimates in the baryon eigenvalues. So we present 
below the result of our calculation, using another set of IK model parameters, given by Gershtein and Dzhikiya 
The wave-functions in this case are given by 

\N >= 0M\ 2 S s > ~0.18| 2 5^ > -0.22| 2 5 A f > -0.051| 4 D M > -0.0039| 2 P A >, (107) 

|A >= 0.98| 4 5 s > +0.16| 4 SV > -0.11| 4 L> S > +0.088| 2 L> M > • (108) 
They correspond to energy values of Mjv = 944MeV and Ma = 1234MeV. The weak amplitudes are: 

Ml = -~^A(k) ^2.1159 + 0.0564fc 2 - 0.00084I 4 ) , (109) 

E2 = -4^® (l0.29 + 1.03fc 2 + 0.236^) x 10~ 3 , (110) 

S v = l -=A(k)k 2 (0.1330 + 0.00171fc 2 ) , (111) 

6V15 V ' 



L v = V3A(k)— (-0.00171 + 0.000549fc 2 ) , (112) 



ni 



[M1) A = -^ A (k) (2.1114 + 0.0567I 2 - 0.00084^) , (113) 
(E2) A = -^ A (k) (0.0521fc 2 - 0.236^) x 10~ 3 , (114) 
L A = -A(k) (0.862 + 0.0189fc 2 - 0.000737^) , (115) 



S A = J- — A(k) (0.337 + 0.00480/c 2 - 0.000172A r4N ) . (116) 
V 3 m V / 

This gives us a feeling of the sensitivity of the weak helicity amplitudes to small changes in the IK wave functions. 
The form-factors calculated from this choice is labeled as "Isgur-Karl 2" in Figs. 1-11. 
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We also include the result of a D-state mixing model, suggested by Glashow pq |, and further discussed by VBJ 
16| , Bourdeau and Mukhopadhyay |32] : 



\N >= VT^\ 2 S S > +Vl\ i D M >, (117) 

|A >= y/1 - 3/3\ 4 S s > ~V^\ 4 Ds > +VP\ 2 Dm > ■ (118) 

The purpose here is to get the correct gA by changing 7, and adjust (3 so that the EMR comes out to be the PDG 
recommended value. We find 



and together they give: 



7 = 0.2048 ± 0.0015, = 0.103 ± 0.011, (119) 



g A = 1-257 ± 0.003, EMR = -0.015 ± 0.004. (120) 



The results of this model is included in Table IV and Figs. 1-11. 

A recent experiment |37j shows that EMR = —0.026. This implies, in the D-state mixing model, (3 — 0.073. For 
brevity, we do not give numerical results for the helicity amplitudes for this case. 

E. The (E2) A / (Ml) A ratio 

The ratio E2/M1 in the 7-/V — > A transition has attracted a lot of attention of theorists and experimen- 

talists J3f| alike. In analogy to this quantity, we can also define the ratio 

AEMR = (E2) A I (M1) A , (121) 

for the axial- vector analogue of the EMR for the process (3). In the SU(6) symmetry limit, it is zero identically. 

We give in Fig. 11a plot of the ratio as a function of Q 2 , keeping W fixed at the Delta excitation, as computed from 
the QM, and compare it with the Adler assumption that it is zero identically, an assumption made in the Kitagaki 
et al. jl8[ analysis. It would be interesting to determine this quantity directly from the experiment in the future. 
We recall here that the EMR , for the vector current, is not well-determined in the QM of the IK variety. While 
Davidson et al. [p5| found the value for this ratio at the photon point to be 

EMR (Q 2 = 0, [25]) = -0.0157 ± 0.0072, (122) 

from the pion photoproduction in the Delta region. New experiments J|7j at the Brookhaven LEGS facility have 
indicated this ratio to be at least twice as large. The IK quark model predicts the ratio at the photon point 

EMR (Q 2 = 0; IK Model) = -0.0033. (123) 

The alternative wave-functions of IK model (Isgur-Karl 2), given above (Eq. (105-106)), yield an EMR of -0.0051. 
Thus, the EMR, obtained in the IK quark model, is much smaller in magnitude compared to the phenomenological 
value. There are indications in the latest Skyrmion approach Jis] ] that the EMR at the photon point is indeed much 
larger than the prediction of the IK quark model, close to -0.05. Thus, the experimental determination of the AEMR 
may be less difficult than what the quark model suggests, as the latter may turn out to be a gross underestimate. 
This quantity should be computed in the soliton models. 

While we are on the subject of the EMR and AEMR, we should make some remarks about the issue of the behavior 
of these ratios as functions of Q 2 , a very topical question. In the perturbative QCD approach (pQCD) [Q, these 
ratios must approach unity, as Q 2 — ► 00: 

EMR (Q 2 00) => 1, AEMR (Q 2 00) => 1, (124) 

while in the IK model, these behave as follows: 

6.26 + 0.627I 2 + 0.252I 4 in _ 3 
EMR = =5 ^ x 10 , (125) 

2.0237 + 0.0623fc - 0.00133/c 
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„ 0.325/c 2 - 0.251I 4 „ o . . 

AEM.R = =5 x 10" 3 . (126) 

2.0204 + 0.0625/c - 0.00133fc 

The last two equations are, of course, meaningless at large Q 2 (that is, large k ). We are giving the last two equations 
merely to indicate very different limits for the EMR and AEMR, at large Q 2 , in the IK quark model, compared 
with pQCD. Nevertheless, it is amusing to note that the IK quark model reproduces the signs and rough orders of 
magnitudes of the pQCD limits mentioned above. 

Several questions are important here: (1) Do the EMR and AEMR approach unity for large Q 2 attainable in 
the laboratory? If the answer is yes, at what Q 2 1 (2) Where does the IK-type quark model predictions definitely 
break down? (3) What happens to the Bloom-Gilman duality jll| in the axial vector sector, on which we have no 
experimental information? Future research should focus on these issues. We may add here that there is a strong 
debate in the literature on the question (1). The pQCD "believers" |i(J tend to see the asymptopia at hand at 
Q 2 > 6 GeV 2 , while the "non-believers" [|42| argue that the Q 2 value for the pQCD rules to be valid is too large for 
us to worry, at the "modest" Q 2 currently available. 



F. Further comments on the comparison of our work with others 



Previous to our work, many authors have investigated the neutrino excitation of the Delta resonance off nucleons, 
but all in the context of the SU(6) symmetric quark model • Schreiner and von Hippel [ [Tcfl have reviewed most 
of these attempts. Among the later works, mention may be made of the work of Abdullah and Close j8|. In this work, 
the CVC condition is imposed, also, quark structure functions are introduced, and the quark form factors are treated 
as constants. In the work by Andreadis et al. again in the SU(6) limit, quark form factors are introduced. We do 
not use form factors at the quark level for the following reason: (1) we want to test quark model predictions rather 
than fitting a constituent-quark form-factor. (2) We believe that the effects of the gluons and sea-quarks have already 
been absorbed in the effective potential of Isgur-Karl quark model. (3) Our results without the form-factor is very 
good compared to the currently available experiments. Thus, there is no phenomenological reason to introduce quark 
form factors. (4) Introduction of quark- level form factors would complicate even more the relationship of the quark 
model to QCD. As it is, this relationship is far from clear. 

Finally, we should discuss the difference between our approach here and a recent work by Hcmmcrt, Holstcin and 
Mukhopadhyay(HHM) |L4| dealing with the NN and NA couplings in the quark model (see Table IV). In that work, 
relativistic corrections to the nucleon gA are taken into account and an agreement with the experimental gA for the 
nucleon is reached. But it yields an off-diagonal NA axial coupling substantially lower than the experimental values. 
In contrast to HHM, the non-relativistic IK approach, used here, does not have the relativistic correction taken into 
account. Thus, the diagonal value of gA still remains off the experimental value (gA(IK) = 1.63). But the off-diagonal 
axial vector matrix elements come out better. Overall, the color hyperfine interaction does not remove the discrepancy 
between quark model estimates and experiments in either approach, when both diagonal and off-diagonal effects are 
computed. Thus, the problem of the quark model in simultaneously explaining diagonal and off-diagonal observables 
does not disappear. 



VII. SUMMARY AND CONCLUSIONS 



We have computed weak amplitudes of the N — > A transition in the framework of the IK quark model, thereby 
dealing with the effect of the color-hyperfine interactions in these amplitudes. Our main conclusions are: 

1. The deficit of the nucleon to Delta (1232) ma gne tic dipole amplitude in the IK quark model estimate, when 
compared with pion photoproduction analysis p5j, is confirmed via the vector form factor at q 2 — > 0. This 
deficit seems to heal around of Q 2 = 0.25 GeV . 

2. There is a mild violation of the CVC in the IK model. Thus, the amplitude , which should be zero by CVC, 
is predicted in the IK quark model to have a small but non-zero value. 

3. The axial-vector transverse amplitudes are largely well-described in the IK model. An exception is our inability 
to get the PC AC value of the Cf. This is not surprising, since we do not have explicit meson degrees of freedom. 
The IK model also violates axial current conservation in the chiral limit. Despite this shortcoming, we are able 
to reproduce the off-diagonal Goldberger-Treiman estimate of the Qf(0) in the IK quark model. 
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4. There is a simple way of parameterizing the SU(6) breaking effects, through a relation (Eq.(56)) that connects 
the two transverse helicity amplitudes to the Adler form factor C^. The IK model gives an estimate of this small 
effect. Its experimental verification, though indirect and difficult, should constitute an important experimental 
challenge, analogous to the determination of E2/M1 in the vector (electromagnetic) sector. 

Hopefully, new weak interactions studies in the nucleon resonance region (1-2 GeV of W) will be possible at 
existing neutrino facilities. There are also new experimental possibilities at the CEBAF on the weak charged and 
neutral current explorations [43-47] of the isobar physics. Though these are intrinsically very difficult, there is some 
hope that such experiments would be possible. They would go a long way towards our understanding of the axial 
vector response of the nucleon, in particular. The role of the Delta isobar, already important in the Adler- Weissberger 
sum-rule j48| , would be interesting to be explored further in the nucleon to Delta weak excitation domain. 
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< A*Ss\N' z S s > 

< A 4 S s \N 2 S s i > 

< A 4 S S ,\N 2 S S > 

< A 4 S S \N 2 S M > 

< A 4 S S \N 4 D M > 

< A 4 D S \N 2 S S > 

< A 2 D M \N 2 S S > 

< A 4 S S ,\N 2 S S , > 

< A 4 S 3 ,\N 2 S M > 

< A 4 S S ,\N 4 D M > 

< A 4 D a \N 2 S s/ > 

< A 4 D S \N 2 S M > 

< A 4 D S \N 4 Dm > 

< A 2 D M \N 2 S 3 , > 

< A 2 D M \N 2 S M > 

< A 2 D M \N 4 D M > 
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TABLE I. Axial-vector longitudinal and scalar amplitudes for various nucleon to Delta SU(6) configurations. Common factor 
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< A*S S \N'S S > 

< A 4 S S \N 2 S S , > 

< A 4 S s ,\N 2 S a > 

< A 4 S 3 \N 2 Sm > 

< A 4 S S \N 4 D M > 

< A 4 D s \N 2 S a > 

< A 2 Dm\N 2 S s > 

< A 4 S S ,\N 2 S S , > 

< A 4 S S ,\N 2 S M > 
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< A 2 D M \N 2 S S , > 

< A 2 D M \N 2 S M > 

< A 2 D M \N 4 D M > 
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TABLE II. Transverse axial-vector amplitudes for various nucleon to Delta SU(6) configurations. Common factor is 
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TABLE III. Vector amplitudes amplitudes for various nucleon to Delta SU(6) configurations. Common factor is 



Y 4 2K^ e h / ' f r ^ ne t ransverse contribution listed is from H X only. Transverse contribution from HY equals to its axial 
counterpart (Table 2) times fc/2m. Contributions from the other configurations not shown are zero. 
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TABLE IV. A comparison of the Adler form factors in the N — + A weak transition in different approaches indicated in the 
text. The momentum transfer squared is taken to be zero. The last four rows are from this work. 
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The Figure Captions 

Figure 1: The ratio A\, % I^A\, 2 vs. Q 2 . 
Figure 2: The ratio A A /2 /V3A A /2 vs. Q 2 . 

Figure 3: C% as a function of Q 2 . Results from the two versions of the IK quark-model, SU(6) limit and the D-wave 
mixing model, are compared with the experimetal results of Kitagaki et al |lS| ]. The experimental errors shown 
here are due to the uncertainty in My (Eq.(93)). 

Figure 4: C\ as a function of Q 2 . 

Figure 5: vs. Q 2 . This form-factor is assumed to be zero in experimental fits. Its non-zero value here indicates 
the violation of the magnetic dipole dominance. 

Figure 6: C 6 y vs. Q 2 . Its non-zero value indicates the degree of CVC violation in the quark model. 

Figure 7: C£ vs. Q 2 . Its non-zero value indicates the violation of magnetic dipole dominance. 

Figure 8: vs. Q 2 . 

Figure 9: C A vs. Q 2 . Results from two versions of the IK quark-model, SU(6) limit and D-wave mixing model 
are compared with experimental results of Kitagaki et al. The experimental errors shown are generated by the 
uncertainty in the fitted parameter Ma (Eq.(98)). 

Figure 10: C§ vs. Q 2 . Since the quark model does not have pion-pole term built in it, this is only the non-pole 
contribution. 

Figure 11: E2/M1 and (E2) A /(M1) A , EMR and AEMR respectively, vs. Q 2 . 
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